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Abstract 

We propose a perturbative-variational approach to interacting fermion systems on ID and 2D 
lattices at half-filling. We address relevant issues such as the existence of Long Range Order, quantum 
phase transitions and the evaluation of ground state energy. In ID our method is capable of unveiling 
the existence of a critical point in the coupling constant at (t/U) c = 0.7483 as in fact occurs in the 
exact solution at a value of 0.5. In our approach this phase transition is described as an example 
of Bifurcation Phenomena in the variational computation of the ground state energy. In 2D the 
van Hove singularity plays an essential role in changing the asymptotic behaviour of the system for 
large values of t/U. In particular, the staggered magnetization for large t/U does not display the 
Hartree-Fock law (t/U)e~ 27T V /t/>u but instead we find the law (t/U)e~^ t t u . Moreover, the system 
does not exhibit bifurcation phenomena and thus we do not find a critical point separating a CDW 
state from a fermion "liquid" state. 
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1 Introduction 



The study of strongly correlated fermion systems apart from its theoretical interest has important 
practical applications ranging from high-T c superconductivity to quasi-one dimensional systems. It 
is usually emphasized the difficulties in handling models describing correlated electrons, such as the 
Hubbard model. This fact has led to various types of approximations schemes, which can be roughly 
classified as perturbative, variational, renormalization group techniques and numerical. The difficulty 
of the task suggests perhaps that one should combine various techniques in order to come as close as 
possible to the exact solution, which is probably imposible to unveil except in one dimension where 
there are models which happen to be integrable. 

The aim of this paper is to combine perturbative and variational techniques in the study of strongly 
correlated fermions. We want to apply to these systems the ideas developed in reference Q which were 
applied to the problem of the Ising model in a transverse field. 

Let us suppose that we are given a hamiltonian of the form H(X) = Hq + XHi, where Hq has 
a nondegenerate ground state ipo and A is a coupling constant. In reference [|lj it was proposed to 
construct the ground state ip(\) of H(\) as 

oo 

V(A) = exp(£ X n U n ) (1-1) 
n=l 

Solving perturbation theory in A to order v implies the knowledge of the collection of operators 
{U n } n= i t „, tl/ . Each operator U n consist in fact of a sum of "irreducible" operators Vi, 

U n = Y,Pn,I V I (1-2) 
I 

Hence inserting ( |1.2| ) into fll.ip and interchanging the order of the sums one arrives to a "dual" de- 
scription of the ground state 

V>(A) = exp(5>i(A)Fi) Vo (1.3) 
i 

where a/ (A) = J2 n ^ n Pn,I- 

This expression suggests an alternative approximation to the ground state ip(\) which consists in 
choosing only a class of irreducible operators Vj whose weights oti are determined variationally. This 
was precisely the approach applied in |Q] to the Ising model in a transverse field, and it is our purpose to 
apply it in this paper to a system of spinless fermions defined on a hypercubic lattice in D dimensions. 
The hamiltonian for this system is given by: 

H = -tj2 (cfe + 4 c *) + E u ( n i - \)( n i - \) (1-4) 

<i,j> <i,j> 

where Cj(cJ) are fermion annihilation (creation) operators satisfying standard canonical anticommuta- 
tion relations (CAR algebra), n« = cjcj is the number operator, t is the hopping parameter and U > is 
a repulsive coupling constant mimicking the residual Coulomb interaction. The sum in ( |1.4| ) is extended 
over all the links of the lattice. The lattice we shall be dealing with is made up of L sites in each spatial 
direction. This amounts to a total number of sites N = L D and it is taken to be even for simplicity. 
The hopping parameter is defined positive t > and Periodic Boundary Conditions are assumed in 
each direction. 

We shall be concerned with the study of the ground state of H when the system is at half-filling, 
i.e., when the number of electrons N e is half the number of lattice sites N e =< J2i n i > = N/2 (so 
that by expanding Ql.4| ) it means that the chemical potential is fixed by [i = DU), and secondly when 



1 



the system is in the strong coupling r egime, i.e., U > t. This means in particular that our unperturbed 
hamiltonian is 



H = ]T U(m - l -){n 3 - l -) (1.5) 

<ij> 

For U > and at half-filling there are two ground states of Hq corresponding to the unperturbed charge 
density states \CDW(e/o)), 

\CDW{e/o) >= J[ cj]0 > (1.6) 

i6A e/o 

where we have taken into account that the hypercubic lattice with N even can be divided into two 
interpenetrating sublattices A e "even" and A D "odd" (or black and blank sites). 

Even though Hq has two different ground states we shall apply our technique to each of the CDW 
states ( |1.6[) . This eventually means that after switching on the hopping term in ( |1.4| ) each initial CDW 
evolve into different CDW-states ip e / a , characterized by different values of the staggered magnetization 
m st , 

™ s * = § < Et- 1 )^ > ( L? ) 

i 

where (— l) 1 = 1(— 1) for i £ A e (A G ). For the initial CDW states take on the values m s t(even / odd) = 

(+V-i). 

An interesting question arises as to whether m s t vanishes for a critical value of (t/U) c different from 
infinity. Equivalently this amounts to the existence or not of Long Range Order (LRO) in the fermion 
system at zero temperature. It is clear that if t/U = oo the ground state of ( |1.4|) is a Fermi-sea of free 
particles which has no preferred distribution in the space which in turn implies vanishing staggered 
magnetization m s t = 0. In 1-dimension the model (|l,4| ) can be solved exactly || and one obtains as 
the critical point the value (t/U) c = 1/2. As a matter of fact, the hamiltonian ( |1.4| ) is equivalent to 
the anisotropic XAZ-Heisenberg hamiltonian, namely, 

1 N 

Hxxz = —„ + ° v y i+1 + Aafa! +1 ) (1.8) 

z %=\ 

This mapping is carried out by the standard Wigner-Jordan transformation and thus the relation 
between the anisotropy parameter A and those appearing in (|0|) is given by, 

The value A = —1 (i.e. (t/U) c = 1/2) corresponds precisely to the Isotropic Antiferromagnetic Heisen- 
berg model in (|1.8| ). The unperturbed \CDW(e/o)) states of the fermion model are mapped onto Neel 
states in the Heisenberg model. We must recall here that mean field theories or Hartree-Fock approxi- 
mations to (|0|) always yield a non-vanishing staggered magnetization m st no matter how small is the 
coupling constant U or the dimension of the lattice. This result is clearly wrong as we have already 
pointed out above for D = 1 and is explained on the basis that the HF method does not take properly 
into account the quantum fluctuations which are very strong in ID. In reference (2| this issue has 
been analized using RG-techniques showing that the existence of a solution of the gap equation for the 
CDW order parameter can be traced back to the fact that in perturbation theory one is disregarding 
some diagrams (those of BCS-type according to the terminology of reference [Q]), and when they are 
properly taken into account they balance the contributions from the diagrams which favor the CDW- 
instabilities. The outcome of this one-loop RG-calculation in D = 1 is a zero beta function (3{U) = 
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for the coupling constant U. This result is by all means only valid in the perturbative regime (i.e. 
U/t small) and implies the existence of a non-trivial fixed point of the Renormalization Group which 
correspond to a Luttinger liquid. Application of the same technique to the case of higher dimensions 
D > 1 always encounters a CDW-instability. 
In summary we have the following situations, 

• (t/U) c = 1/2 in D = 1. Exact. 

• (t/U) c = oo in Mean Field Theory, MD. 

• (3{U) = in Renormalization Group Approach, D = 1. 

We shall address in detail the issue of the critical value (t/U) c in the system of interacting fermions as 
well as other problems that emerge in connection to this one. 



2 The Perturbative- Variational Ansatz for Interacting Fermions 

According to the perturbative-variational (PV) method developed in reference (l| we must first deter- 
mine the form of the set operators {U n } by inserting the exponential ansatz ( |1.1| ) in the Schrodinger 
equation for H. Then these operators serve us to construct variational wave functions by inserting 
them back in the exponential ansatz. The perturbative equations that determine the U n for the lowest 
order in perturbation theory obey the equations, 

order A: ({H , U x ] + ffi)Vo = £ (1) V>o (2.10) 



order A 2 : ([H , U 2 ] + [H u CT X ] + 1/2([[H , U X \M))^ = £ (2 Vo (2.11) 
where and E^ are the perturbative energies to first and second order. 

A solution of equations ( [2.10 ) and ( 2.11|) in the case where Hq is given by equation (^) and H\ is the 



hopping term of ( |1.4| ) is the following, 



Ui 



2D 

<h3> 



I E T a ( 2 - 12 ) 



U-2, — E Cijk{TijT k i (2-13) 

<i,j;k,l> 

where Tij = c\cj + CjCi and < k, I > is a path on the lattice from the point i to the point / without 
* 7^ j 7^ k 7^ ' an d Cijki some constants which may depend on the shape of the path and on the 
dimension D, whose precise value is not important for our discussion. The perturbative energy up to 
second order is 

NUD, 4 . t , 

E = 1 + (-) 2 (2.14) 

4 v 2D - VU J ' v ' 

Observe that U\ is proportional to the kinetic term in the original hamiltonian H. The solution of 
these equations is by no means unique. We can add to XJi an operator of the form Tu with i, I on the 
same sublattice. The action of such an operator on the state \CDW) is zero. 

Now the second stage of the method, according to eqs. (|l.l| ) - (|l.3|) , consists in exponentiating the 
U n operators ( 2.12j ) and ( |2.13| ). The simplest ansatz takes into account only the U\ operator yielding 



V$(a) = exp(| E {c\ C] + t%*))\CDW (e/o) > (2.15) 

<ij> 
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whereas if we take into account the second order operator U2 the ansatz will become 



^(a) = exp(| J2 T ij + \ E *iW T i3 T kl)\CDW (e/o) > (2.16) 

<i,j> <ij;kl> 

where ctijki will only depend on the shape of a 4-step non-backtracking path. In both cases the 
parameters a's are used to minimize the energy of the trial states, which fixes their dependence on the 
ratio t/U. For t/U small we must have a ~ t/U, ocijhi ~ (t/U) 2 and such that (ip^) agrees to 
first (second) order with the exact ground state of H. 

We notice that this is the essence of the PV method proposed in Q, namely, to use the results obtained 
to a given order in perturbation theory in order to establish a variational calculation which extends the 
validity of the perturbative results beyond perturbation theory. 

In what follows we shall restrict ourselves to the analysis of the first order ansatz ip^^ct) (hereafter 
called tp(a)). 

As t/U increases we expect also a to increase, after all the role of this hopping term in ( |1.4| ) is to 
disorder the CDW state, but that is precisely what the exponential operator in the ansatz fl2.15|) is 
doing on the CDW-states. 

Eventually, as t/U increases we may get a = 00 in which case, as is shown below, the state ( p. 15 ) 
becomes the Fermi-sea or ground state of H\, 

lim ip e / (a) ~ iFermi Sea >= IT ct|0 > (2.17) 

a— >±oo ' 

kGFermi Sea 

All the states ip{ a ) with a finite have non- vanishing staggered magnetization and hence only for a = 00 
we get a state with m s t = 0. This in turn implies that the value of t/U at which a = 00 is nothing else 
but the critical value (t/U) c . Therefore the rule to find in our approach the critical value (t/U) c is given 
by finding the minimum of the variational energy E(a;t/U) of the state ifj(a) for t/U = const.. This 
fixes t/U as a function of a, i.e., t/U = !F{a). Then the critical value (t/U) c at which the staggered 
magnetization vanish is given by 

(t/U) c = lim T{a) (2.18) 



If (t/U) c turns out to be non-zero then the ansatz ( 2.15| ) would be valid in the whole range < t/U < 



(t/U) c and beyond this the variational state remains the free Fermi sea ( p. 17 ). Hence to get more 



information we should study the model from the "other" side of coupling constant, as for example 
perturbation theory in U or some other method. 

To prove ( 2.17 ) notice that the ground state of the operator appearing in the exponential of ( [2.15 ) 



is, for either sign of a, given by the Fermi sea at half-filling. Now if we let a go to 00, and we introduce 
the resolution of the identity we project the state ip(a) onto \Fermi Sea). We note that a similar 
technique is used in the so called projected Monte Carlo method in order to filter out the ground state, 



lim e~ tH \<t> >~ \G.S. of H >< G.S. of H\<j> > (2.19) 



t— »CXD 



It is assumed that the trial state \<j) > is not orthogonal to the true ground state of H. 

For the ansatz ( p. 15 ) we do not have to resort to numerical computations to obtain the mean value of 



H since this state is in fact a Hartree-Fock state of a restricted type. 



To prove this we shall express the state ( 2.15 ) in momentum space. For that purpose we introduce the 



creation operator c\ and annihilation operator Ck which are the Fourier transformed of Cj, 



1 ^e ih - R >4 (2.20) 
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where k = x( n i' • • • > n D), with rtj = 0, 1, . . . , L — 1 mod(L). 
The CDW states (|1.6| ) can then be expressed using Q2.20|) as, 

\CDW(e/o) >=i=n(4±4 +Q )|0> 



(2.21) 



where Q = (n, . . . , 7r) is the AF vector and ]J k stands for the product over the reduced Brillouin zone, 
also called magnetic zone. 
Taking into account that: 



J2 ( c i c i+ c j c i) = Y, e k n k 

<i,j> k 



D 



-2^ COs(£; a ) 
a=l 



(2.22) 
(2.23) 



and rik = c\,Ck, we can write the states ( |2.15| ) as 



V'e/oH = -^ne-f £ «=4±e-W«4 +Q )|0 



> 



(2.24) 



From this expression we can inmediately derive equation ( |2.17 ). 
Concerning this result we would like to make some comments: 



i) Equation ( [2,24 ) shows that the PV ansatz ( [2, 15 ) is a Hartree-Fock state, i.e., it is a state of the form 

(2.25) 



ipHF = Y[(u k cl + v k cl +Q )\0 > 
k 



with 



(2.26) 



ii) Since the functional dependence of u k and v k is fixed by equation (2.26) he PV state is a restricted 
Hartree-Fock state. This implies that the variational method applied to the PV state using a as 
the unique variational parameter, needs not coincide with the variational method applied to an 
unrestricted HF state where all the parameters Vk are variational Q, |7j. 



The computation of the mean value of the hamiltonian H (]1 .4f ) on the state ( 2.24j ) is standard, we 
quote here the result and leave the details of the proof for the appendix A, 



where e(t/U,a) is the reduced energy 



E = —-NDU e(t/U,a) 



t 



e(t/U, a) = 4~ii(a) + / 2 (a) + if (a) 



(2.27) 



(2.28) 



and Io(a), h(&) are the following integrals, 
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. . f d D k 1 

M« = / , 9 vd u \ 2.29 

Jn (2tt) 1j cosh(ae fc ) 

1 f d D k 

= 2D In J^D e k tan H^k) (2.30) 



Here Q denotes the interior of the Brilluoin zone (— 7r < A: a < tt). 
Minimization of E with respect to a yields 



t/U = -\{h + I -ft) (2.31) 



where I = dio/da, etc. For a given value of t/U there may exist several values of a satisfying ( p. 31 ) 



The desired value is selected by minimizing the energy E. This equation will play the role for us of the 
standard gap equation of the HF method. 

We can perform a consistency check of our formalism which consists in verifying that for small 
parameter a it is linearly related to the effective coupling constant t/U. This must be so because 
it means that we are very close to the ground state of Hq, the CDW state, and thus we are in the 
perturbative regime where both quantities are proportional, for that is the origin of the construction 
of the U\ operator (|2.12j ). 



To prove this assertion we must evaluate the leading behaviour when a — > of the basic integrals Iq, 
1 1 in ( 2.29j ), ( |2.30 ). From these expressions we arrive at, 



Iq ~ 1 Iq ~ —IDOL 



T a^O T , a— >0 .. 

1\ ~ a i 1 ~ 1 



(2.32) 



Inserting these small-a behaviours in the defining relation ( |2.31| ) of t/U in terms of a we find 

t/U ~ a(D - -) (2.33) 

which is indeed the linear relation previously advanced in equation ( |2.12j). This relationship is made 
explicit in Figure 2 where we have plotted the corresponding equation ( |2.33| ) for the cases D = 1 and 
2 dimensions. 

Once this check is done the big issue is to see whether we can extract nonperturbative results for the 
system of interacting fermions by going over large values of a for which the perturbative relation ( [2.33 ) 
ceases to be true. In the forthcoming sections we investigate this issue in one and two dimensions. 



3 Strongly Interacting Fermions in D = 1 and the Heisenberg Model 

Let us consider the system of spinless fermions interacting by means of the hamiltonian H (|1.4j ) in 
one dimensions. This is a very special case where exact formulas exist for the magnitudes of interest: 
ground state energy, magnetization or Long Range Order, density correlators, etc .... This allows us to 
confront our PV methods predictions with the exact results. Moreover, this is a very demanding test 
ground for it is well-known that mean field theory treatments are usually doomed to failure because 
the strong quantum fluctuations occuring in one dimensions are not properly taken into account. 

In one dimension the system undergoes a quantum phase transition at the critical value of the 
coupling constant (t/U) c = 1/2. This is an essential singularity phase transition. For small values of 
t/U the system is in an ordered phase in coordinate space corresponding to a perturbation of the CDW 
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state. In this phase the staggered magnetization calculated by Baxter is non-vanishing. For large 
values of t/U the system is in a disordered phase similar to a fermion "liquid" of Luttinger type. 

Elucidating the existence of a finite critical value (t/U) c in this case with our method goes over the 
determination of the asymptotic behavior of integrals Iq, I\ and its derivatives for a — > oo. 

Let us outline how to compute the leading large-a behaviour of Iq D ~^ for it reveals the role played 
by the Fermi points in the physics of the system, and it is also useful to generalize the same technique 
to the more involved case of two dimension discussed in next section. 

For a — > oo only the regions in the neighborhood of the points satisfying cos k = 0, that is the Fermi 
points k = it/2 and — vr/2, contribute significantly to the integral Iq in equation ( 2.29| ) for D = 1. 



As matter of fact, both Fermi points contribute the same amount to the integral under consideration. 
Thus, it is convenient to introduce a cut-off A to isolate the contribution around the point, say k = ir/2: 



4 D=1) ~ 2 



f + A dk 1 



A 27r cosh(2a cos k) 

Shifting variables k = ir/2 + k with k small and thereby approximating cos k ~ — «, we may write 

2 f A , 1 



jf =1) „ ^ / dK 



2tt J-a cosh(2a/-c) 

Changing variables again x = 2an we can freely take the limit a — > oo in the limits of integrations and 
arrive at 

r (2?=i) 1 f°° dx 1 

In ~ / 

2ira J-oo cosh x 2a 

This result establishes that Iq D ~^ vanishes as a power of a with exponent —1 in the asymptotic region 
a — > oo. It is straightforward to extend this analysis to the other integrals I , I\ and I[ so that we 
summarize the results hereby: 

(3.34) 

1 7T 1 24t? 



Now, in order to determine the critical value (t/U) c as prescribed in equation ( 2.18| ) we must insert 



the leading values ( |3.34 ) in the defining variational relation of t/U ( |2.31| ) This inmediately yields the 



critical value (t/U) c = - ~ 0.63662 which is quite close to the exact value (t/U) c = 0.5. 

However this quick derivation is not correct. The first indication that something peculiar is going on 



in one dimensions comes from the shape of t/U as a function of a in equation (2.31). This function 
has the horizontal asintota 2/tt as explained above but it is not a monotonous increasing function. In 
fact, it exhibits a local maximum at a = 2.3307. This in turn means that the relationship ( |2,31| ) is not 
invertible so that we cannot express a as a function of t/U to be inserted in the energy equation and 
the remaining formulas. 

Therefore, to determine the correct critical value (t/U) c we must proceed more carefully and in doing so 
we shall unveil the peculiarities of the PV method when applied to the special case of one dimensions. 
The correct procedure is exemplified in Figure 1 where we plot the reduced energy e^ D=1 \a; t/U) as a 
function of a for varying values of the coupling constant t/U which is considered as a parameter. In this 
fashion we are correctly implementing the minimization program thereby assuring that we are picking 
up the absolute maximum of e^ D=1 \a;t/U) (recall that a maximum in the reduced energy correspond 



to a minimum of the ground state energy E ( p. 27 ) ), not just the local maximum. This local maximum 



contribution was the one we were only considering previously when the critical value 2/tt was naively 
obtained. 
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For small values of t/U, say 0.2 (see Figure la), the function e^ D=l \a;t/U) exhibits an absolute 
maximum for a finite value of a and an absolute minimum for a —* oo. 

When t/U is increased, say 0.4 (see Figure lb), the relative distance between the maximum and the 
minimum is reduced. 

Eventually, for say t/U ~ 0.735 (see Figure lc), a local minimum shows up nearby the absolute 
maximum and the minimum at a —* oo becomes a local maximum. Yet, the first maximum dominates 
over the new maximum at infinity and it is thereby picked up in our minimization procedure. 
We are able to compute the exact value taken by e^ D=1 \a; t/U) in the oo, as a function of t/U. In 
fact, substituting the asymptotic behaviours (|3.34 ) in e^ D=l \a;t/U) we get, 



e^\oo;t/U) = ^ + {-f (3.35) 

7T U 7T 

This is an increasing function of t/U and eventually this local maximum will overcome the thus far 
absolute maximum occurring at a finite a. When this fact happens, we are obtaining the true value of 
the critical point (t/U) c . 

As the coupling constant t/U keeps being raised, two things happens: the local minimum moves towards 
the yet absolute maximum and the relative height of the maximums is reduced. The critical point occurs 
when the height of the local maximum at a infinity equals that of the maximum occurring at finite a. 
This happens for the following value 



(t/U) c = 0.748375 (a c = 1.80032) (3.36) 

From this value of t/U the maximum at infinity becomes the absolute maximum henceforth. 
It is worth pointing out that in our procedure for one dimensions explained above we have been implicitly 
describing an example of what has come to be known Bifurcation Phenomena in several branches of 
Physics (We shall see in next section that such phenomenon is absent in two dimensions.) 
As a matter of fact, this way of thinking is well suited for our problem as well as enlightning. Bifurcation 
phenomena typically occurs in a system subject to the action of a certain potential say V(x;p) where x 
can be a generalized coordinate describing the state of the system and p an external parameter which is 
somehow varied. For a certain initial value of p, the potential V exhibits two minima x m \ < x m 2 with 
V(x m i) < V(x m 2) and the system locates at x m \. As the external parameter p is varied, the two minima 
start approaching each other till eventually they meet, V(x m i) = V(x m 2). This in turn defines the 
critical value p c . As p is yet more increased, the two minima exchange their roles, V(x m i) > V(x m 2), 
and the system locates at x m 2- Altogether, when p is varied at one go over its full range of validity, the 
system undergoes a "phase transition" from the state at x m \ to the state at x m 2- But this is precisely 
what we have been describing above in our method where e^ D=1 \a; t/U) plays the role of the potential 
function V(x;p), a ~ x and t/U ~ p. Interestingly enough, the phase transition we are describing in 
this fashion is the quantum phase transition occurring in the spinless interacting fermion system from 
the CDW state to the Fermi sea state. 

This result is quite different from what one obtains using mean field theory, where there is always a 
solution of the gap equation for Acdw n ° matter how small is the coupling constant U. The physical 
explanation of this result is that any perturbation having momentum tt links the ground state, given 
by the Fermi sea at half-filling, to states of arbitrarily low energy in which a particle just below right 
(left) Fermi point is pushed above the left (right) Fermi point. In this way the Fermi sea at half-filling 
becomes unstable under this perturbation turning into a CDW state. Our result above in D = 1 is 
precisely the opposite. It shows that starting from the CDW-ground-state it becomes unstable for t/U 
bigger than a critical value (t/U) c . 
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We can summarize our discussion in the following diagram, 



Mean Field VC/,VL> 
Fermi Sea — ► CDW 

PV-Method D = 1, t/U > (t/U) c 
CDW — ► Fermi Sea 

Of course the correct result is that the CDW state for t/U > 1/2 should become a Luttinger liquid, 
but it is clear that in our approach we have not introduced enough parameters to tell the difference 
between a free Fermi liquid and a Luttinger liquid. 

The outcome of the procedure oulined in Figure 1 is to achieve the correct set of values (t/U, a) 
associated to the PV method. When this collection of values is plotted it looks like in Figure 2. 
Notice firstly that t/U as a function of a is single- valued and secondly, that from a finite initial value 
a c = 1.80032 the critical value (t/U) c = 0.748375 is reached once and for all. 

Let us now turn our attention to physical quantities and consider first the staggered magnetization 
m s t fll.Tp or equivalently the averaged fermion density depending on whether we choose the anisotropic 
Heisenberg viewpoint (1.8) or the spinless fermion description ( |1.4| ). 



First of all, notice that the staggered magnetization m s t (|1.7| ) precisely coincides in our method with 
the basic integral Iq(ol) ( [2.29 ) 



m st (t/U) = I (t/U) (3.37) 
In fact, by going to momentum space and using the results of Appendix A we find 

2 2 1 

mst = N < * >= N £ < C *° k+Q >= £ c^o^) 

However, in order to extract the physics underlying this magnitude we must not plot Iq as a function 
of a but instead as a function of the coupling constant t/U. After all, the parameter a is somehow an 
artifact of our method which is related to the physical parameter t/U by the plot in Figure 2. Expressing 
a as a function of t/U in Io(a) we obtain the plot of the magnetization m s t which we show in Figure 
3. In this figure we also plot the exact magnetization as given by Baxter [1C]. Several remarks are 



in order. Only for small values of t/U there is a good quantitative agreement between both methods. 
Beyond these values t/U ~ 0.2 there is major quantitative disagreemet. Yet, our approximate method 
gives a reasonable qualitative behaviour in the sense that it predicts the existence of a critical value 
(t/U) c beyond which the magnetization totally vanishes. However the PV method (to lowest order) 
predict an otherwise absent finite jump in the magnetization at the critical value. 

In all, the physical meaning of the magnetization graphics is that we can predict with our method 
the absence of LRO for small t/U at zero temperature for the system of spinless interacting fermions. 

Now we move on to consider another relevant physical quantity such as the ground state energy 
( 2.27] ) as a function of t/U. In Figure 4 we plot the reduced energy ( |2.28 ) as obtained by the PV method 



against the exact result given by several authors [11], [12|. The agreement is considerably much better 



than with the magnetization. In fact there is a very good quantitative agreement not only for small 
t/U but for the whole physical range till t/U = 0.5. The major difference though is the prediction of a 
bigger (t/U) c than the exact result. 

As a precise quantitative check of our values we have collected our PV result for the reduced energy in 
the Isotropic Heisenberg Antiferromagnet (t/U = 1/2) in Table 1 along with the exact result |]] and 
several other methods. The result is a good estimation of the ground state energy. 
Furthermore, we are also able to compute exactly within our approximate method any multi-density 
correlators. For instance, the density-density correlator < (rii — l/2)(rij + i — 1/2) > computed in 
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Appendix A has the interpretation of the Short Range Order in the Heisenberg model for it is equivalent 
to the spin-spin correlator < SfS? + i >. In Appendix A we show how to compute it as a function of 
t/U. We do not show the plot of this quantity against the exact value but only remark that there is a 
good quantitative and qualitative agreement between them. 

Finally, we close the analysis of the D = 1 dimension case by pointing out that the reduced energy 
e(t/U) exhibits a finite jump in its derivative at the critical value (t/U) c = 0.748375. The constant 
value of the derivative at the right hand side of the critical point comes from the linear contribution in 
equation ( |3.35| ) of the absolute maximum at a — > oo previously discussed. The meaning of this fact is 
that our PV method predicts a first order quantum phase transition for the fermion system. 
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4 Strongly Interacting Fermions in D = 2 Dimensions 



After studying the pros and cons of our PV approach to one-dimensional interacting fermions, we 
address the much more unknown case of D = 2 spinless fermions. This case is particularly interesting 
in physical applications such as high-T c superconductivity. 

First of all, we should quote a major difference ocurring with the D = 2 fermions: there is no 
connection with the Heisenberg model in 2 dimensions (as far as the present knowledge is concerned). 



That is to say, the Wigner- Jordan mapping transforming ( |1.4| ) onto (1.8) ceases to work in higher 
dimensions. Correspondingly, all of our conclusions will be only valid for the fermion system but not 
for the spin system. 

The system of 2D interacting spinless fermions is worse understood in comparison with the 2D 
Anisotropic Heisenberg model. The issue of LRO has been elucidated for the spin system in references 
]l3| , ]14| ] using the reflection property in connection with the Peierls argument. It turns out that in 2 
dimensions there exists LRO for low temperature T and small enough anisotropic parameter. The main 
obstacle in extending similar results to spinless fermions relies on the fact that no type of reflection 
positivity property is known to hold for fermions. Recently, though, what has been possible to prove 
the most is the existence of LRO for low but finite temperature and small coupling constant t/U fllfl . 
There is not exact proof of this statement at zero temperature, although it is believed that LRO persists 
for t/U < 1. In fact, numerical methods such as Quantum Montecarlo fails to provide a definite answer 
at zero temperature [^] . 

With this considerations in mind the study of exact properties for two-dimensional fermions at T = 
remains an open problem and we undertake this study with our PV method developed in previous 
sections. 

We follow a parallel analysis as described in one dimensions. The first key observation is that when 
plotting the function t/U = t/U(a) as given in equation ( 2.3l| ) NO local maximum shows up as is the 



case in D = 1. In fact, Figure 2 shows this dependence which turns out to be monotonically increasing. 
This property guarantees that the function is invertible so that we can trade the a parameter in our 
PV results for the more physical parameter t/U. 

Moreover, this conclusion is further backed by a similar analysis as in Figure 1 for the reduced energy 
e(a;t/U) for increasing values of t/U. We have performed such analysis with the result that there is 
no Bifurcation Phenomena present in D = 2 spinless fermions. In other words, the function e(a;t/U) 
always (\/(t/U)) has a maximum at finite a, which in turn is the absolute maximum for the whole range 
of a. The value e(oo; t/U) at infinity never overcomes that maximum neither a local minimum appears 
in the vicinity of that maximum. 

Thus, the minimization procedure of the PV method is entirely governed by the maximum at finite a 
and the analysis is thereby simplified in D = 2 dimensions. 

The second key observation present in Figure 2 is that there is not a finite critical value {t/U) c of 
the coupling constant towards which a tends asymptotically at infinity. Quite the contrary, t/U always 
grows as a function of a. In fact, we can compute the behaviour of (t/U)(a) in the asymptotic regime 
a —* oo for the D = 2 case. To this purpose, we need to find that behaviour for the integrals Io( Q )> 



I()(a), Ii(a) and I'i(a) and then insert them in ( 2.31 ) 



We have left out the calculation of the leading order behaviour of those basic integrals for the 
Appendix B. We hereby present the final results, 

T {D=2), x 

(4.38) 

In doing this calculation it is essential to deal with the van Hove singularities (saddle point singularities) 
appearing at the corners A(n, 0) and 5(0, tt) of the Fermi Surface (FS) for two-dimensional fermions 



In a 
2-na 


^(«) 


In a 

2-Ka' 1 


(I) 2 


da (« 


\ In a 
' 24a 3 
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as in Figure 5. Analogously as in D = 1 dimensions, only the region of the Brillouin zone near the FS 
contribute significantly to the integrals when a — * oo. Much care is needed to isolate this asymptotic 
behaviour as we have to distinguish the contribution of the van Hove points from the rest of the FS 
points which we call regular points (see Figure 5). In fact, in the intermediate step, there is a highly 
non-trivial cancellation of cut-off singularities between van Hove and regular-points contribution. 
Moreover, we have performed numerical analysis of these integrals in the asymptotic region t/U — > oo 
and we find numerical agreement with the results expressed in equation ( 4.3EQ . 



The main feature of equations (4.38) is the presence of the logarithmic terms which were absent in 



D = 1 dimensions ( |3.34 ). This is ultimately the signal of the van Hove singularity and the reason why 



t/U grows indefinitely: substituting equation ( f4.38| ) in t/U ( |2.31 ) we get 



2 3 

t/U ln « ( 4 -39) 



which explains the logarithmic growth of t/U present in Figure 2. 

Therefore, our method applied to D = 2 spinless fermions does not predict the existence of a quantum 
phase transition from the CDW state where we started out to a new phase, say a disordered fermion 
"liquid" state, as was the case in D = 1 dimensions. 

The analysis of physical properties such as magnetization and ground state energy goes in parallel 
to the D = 1 case of section 3. We show our results in Figures 3 and 4. The main feature of these plots 
is the absence of a critical point (t/U) c which makes the magnetization decrease to zero at infinity and 



the reduced energy to grow to infinity on the contrary. In fact, from (4.39) the asymptotic relation 
between a and t/U is 



a ~ eir' /c/ (4.40) 
Upon substitution of (|4.40|) in the magnetization ( 3.37 ) we obtain, 



m st ~ (t/U)e-^ t/u , t/U -» oo (4.41) 

It is interesting to compare the result ( |4.41| ) with those predicted by a Hartree-Fock method |J17|], which 
in 1,2 and 3 dimensions are given by 

m st ~ Tje- 2n tf (# = 1,3) 

(4.42) 



m st ~ ±e~ 2n Vu (D = 2) 

The square root behaviour in the exponential is due to a combined effect of the nesting of the D = 2 
Fermi surface at half-filling and the van Hove singularity of the density of states. Hence, the HF result 
gives a stronger tendency to a CDW state in D = 2 as compared to D = 1 or 3. On the contrary, 
the PV method when extrapolated to the U weak coupling regime, gives a behaviour similar to that of 

2 

D = 1,3 HF, the only difference being in the factor ^- instead of 2-k. It is an interesting question to 
discriminate which of the two behaviours, namely the square root law or the linear jj law in the 
exponent controls the gap in D = 2 dimensions at half- filling. We want to stress that formula ( [4.41 ) 



is obtained from an extrapolation of the PV state ij)(ot) from the perturbative region of the ansatz 
(t/U <C 1) to the non-perturbative (t/U ^> 1) by means of the variational method. We find indirect 
support for this method when taking the t/U — > oo in the PV formulas for we find that the physics 
of the system is governed by the momenta very close to the Fermi surface, and the Fermi surface is a 
concept which makes sense for a Fermi sea located at t/U — > oo in the coupling constant. 

In reference |l7|| a numerical analysis has been done for the 2D Hubbard model showing that the 
e -2nt/u ^ g no ^. inconsistent with the numerical results. 
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Let us notice that in the above discussion we are comparing asymptotic behaviours concerning the 
staggered magnetization. Another similar analysis could be done as far as the gap Aqdw is concerned. 
In the HF method both quantities are easily related by Acdw = Um s t, but in our PV method that 
relationship is not that easy to obtain for one has to go through the computation of the gap by doing 
an extension of the PV ansatz to the first excited state of the system. We shall not dwell here upon 
this and other extensions which will be reported elsewhere. 



5 Final Considerations 

In section 2 and appendix A we have shown that the trial state iftQ) based on the lowest order PV 
method is equivalent to a Hartree-Fock state of a restricted type. An straightforward generalization of 
the ansatz ( |2.15| ) is given by 

ipUHF = exp(i 4 A ijCj) \CDW > (5.43) 

<i,j> 

where A^ is an hermitian matrix (i.e. A*- = Aji). If the matrix Aij is an hermitian matrix only 

depending on the difference Ri — Rj, then it can be shown that ipUHF is an unrestricted HF state of 
the form ( |2.25| ) where the parameters u& and take generic values depending on the eigenvalues of 
the matrix A. 

Carrying out the PV method to highest orders in the coupling constant t/U we find operators as 



the ones in the exponent of ( 5.43 ), along with operators with a more complicated structure. This fact 



may already be seen from the expression of the operator XJ% ( 2.13j ). First of all, XJi is a sum of operators 



involving 4-fermion terms, which implies that the corresponding state ( 2.16| ) is not a HF-state. 
In fact, working with this type of states is not simple at all. The interesting feature of t/^ 2 ) is that it 
contains correlation effects which are absent in the HF-state (A similar phenomena happens for 
the Ising model in a tranverse field ITF where the analoge of the ip^ state is also a mean field state, 
while the second order ansatz is not mean field and it incorporates correlations ). 
From a more physical point of view the effects of the operators U n or Vi is that when they act on the 
unperturbed ground state, say \CDW(even) >, they change locally the even-CDW vacuum into an 
odd-CDW vacuum. From the form of the Ui operator we see that the irreducible operators Vi can be 
labelled by clusters of dimers, i.e., products of bilinear operators T^j = c\cj + CjCi with the links < i, j > 
belonging to some cluster. The linear clusters are reminescent of the string operators considered in 
references but we should in principle consider all possible clusters in order to construct the 

exact solution via the exponential ansatz. This is probably beyond our present possibilities so that one 
has to restrict the class of operators to build an ameneable ansatz. If this constructions were possible, 
the parameters aj would give us information about the size of the CDW-patches in the true ground 
state. Work in this direction will be reported elsewhere. 

In this paper we have chosen a variational technique in order to compute the parameters aj of the 
exponential ansatz. There is however an alternative way to compute this parameters which is called 
the Fokker-Planck method |20[| . This method tries to generalize to quantum lattice hamiltonians the 
well-known fact that in quantum mechanics, given a wave function of exponential form i[)(x) = e~ s ^ 
one can find a hamiltonian Hpp, called Fokker-Planck hamiltonian, for which ip(x) is an exact ground 
state. 

In the case of the ITF model, the Potts model and others the FP-method can be easily generalized 
for states of the form (|1.3| ) due to the fact that the operators Vj commute among themselves (in a sense 
they play the role of the x-coordinate in the quantum mechanical case.) Then, given the ansatz (1.3) 



one can find the corresponding FP-hamiltonian and can in turn compare with the original hamiltonian. 
From this comparison one finds the a'jS in terms of the coupling constants in H. 
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It is an interesting question to see whether the FP-method can be extended to problems with 
fermions. The first problem that we encounter is that the operators Vi no longer commute among 
themselves (see Appendix C). The advantage of the FP-approach is that one can introduce and do 
calculations with a large number of parameters uj which serve to introduce much of the expected 
physics in the ansatz. 

To summarize, we have employed a Perturbative- Variational (PV) method to study strongly in- 
teracting fermions at half-filling and zero temperature. To lowest order in PV we have stressed the 
similarities and differences with the mean field theory or Hartree-Fock method. To higher order in PV 
they clearly differ. In ID our method is capable of unveiling the existence of a critical point in the 
coupling constant at (t/U) c = 0.7483 as in fact occurs in the exact solution at a value of 0.5. In our 
approach this phase transition is described as an example of Bifurcation Phenomena in the ground state 
energy. In 2D the van Hove singularity plays an essential role in changing the asymptotic behaviour of 
the system for large values of t/U. In particular, the staggered magnetization for large t/U does not 

/ IT 2 

display the Hartree-Fock law (t/U)e~ 27r v t / u but instead we find the law (t/U)e~~ t ^ u . Moreover, the 
system does not exhibit bifurcation phenomena and thus we do not find a critical point separating a 
CDW state from a fermion "liquid" state. 

Many of the ideas and techniques that we propose in this paper can in principle be applied to 
the Hubbard model in the strong coupling regime (U > t), which is believed to be one of the most 
prominents candidate for a microscopic description of the high-T c cuprates. 
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A The PV Formulas for Interacting Fermions 



In section 2 we have briefly explained the fundamentals of the perturbative- variational method and at 
the end of it we have set up the basic equations ( 2.27[ )-( ^31 ) which encode the essence of our method 
as applied to the system of interacting fermions. In this appendix we give the highlights and some 
details of the computations leading to those equations. 

The hamiltonian (1.4) which is the subject of this study is split into two pieces, 

H = H + H 1 (A.l) 

F = U £ K - l -){n 3 - \) (A.2) 

<i,3> 

Hi = -t (4 c j + c ) c i) = t^tknk (A.3) 

<i,j> k 

The norm of the variational state 

^(a) = exp(| ]T (ct Cj + c}«))\0 > (A.4) 

<i,3> 

plays a major role in the calculation of all the averaged quantities and correlators of the theory. We 
shall denote this norm by Z for it recalls the partition function of an associated statistical model. To 
see this, let us compute its value. 

Z =< ip(a)\ijj(a) >=< CDW\exp(-aJ2e k n k )\CDW > (A.5) 

k 

To compute this partition function we note that the variational state (A.4) can be given the following 
representation in momentum space, 



iv> >= IIk4 + «fc4+ol° > ( A - 6 ) 

k 



e 2 €fe e 2 6k +Q e 2 efe 

Uk = ^ Vk = ^^ = ^ (A - 7) 
This representation allows us a readly derivation of the norm in terms of the functions u k , v k namely, 
Z = n'k( u k + v k) = IIa: cosh(aefc), and we have made extensive use of the identity e k = ~^k+Q- I n 
order to perform the continuum limit it is convenient to introduce the free energy /jv of the associated 
statistical model by means of the identity 

Z =< il)(a)\il>{a) >= exp(7V/iv(a)) (A.8) 
Then, from the previous considerations we arrive at the final expression for the free energy. 

exp(Nf N (a)) = exp(i ^ ln(cosh(ae fc ))) (A.9) 
1 k 

If dPk 

f(a) = lim f N (a) = - D lncosh(ae fc ) (A.10) 

Interestingly enough, equation ( A.lOj ) represents the free energy of a system of free fermions placed on 
a D-dimensional lattice as we already advanced. 
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The expectation value of the Hi part of the hamiltonian is easy to evaluate in terms of the asociated 
free energy 

dZ 

— = - < iP( a )\J2e k n k \ip(a) > (A.ll) 

k 

< H 1 >=-t^-(lnZ(a)) (A.12) 
oa 

with the following result already expressed in the continuum limit, 

tN f d D k 

< Hi >= — — ^-p-e fc tanh(ae fc ) (A.13) 

This much for the mean value of Hi in the exponential ansatz. As for the Hq part of the hamiltonian, 
we need the expectation value < J2<i j> n i n j > which represents the correlation function of two density 
operators. Transforming the density operator to momentum space 

"* = Jj E e~^- k ^c\c k2 (A.14) 
fci,fc 2 

the afore mentioned density-density operator reads, 

1 D 13 

H n i n i = Jj E 4 1+fc2ife+fc4 ^e^ 3 4 i Cfc 3 4 2 c fc4 (A.15) 

<i,j> ki,k2,k3,k4 a=l 

where k^ 3 = k a \ — k a ^. At this stage it is aparent that the way to compute this correlator is Wick's 
theorem. Thus the fermion correlator in fc-space is 

< 4i c te >= ( S ki-k2,Q + h-,-k 2 ,Q) (A. 16) 

u k 1 v k 2 

where u k = v k+ n. Therefore the basic contractions are the following two, 

< n k >=< ctcfc >= (A.17) 

< cick+Q >= (A.18) 

The contraction in equation ( |A.17| ) represents the average number of particles with momentum k 
and we notice that it exhibits the standard Fermi-Dirac distribution form with the parameter a playing 
the role of the inverse temperature but recall that all our formalism is at zero temperature. 
A lengthy but straightforward calculation yields 

Tf E < W >= 7 ~ T'o " E -A (A.19) 



where 



N 4 4 u ^4 

<?,J> a=l 



1 1 r d D k 1 

o(«) - Jf 2_, cos h(ae fc ) ~" Jn (2tt) d cosh(ae fc ) ( ' 20) 

Ia{a) = ly cos (y (A.21) 

k 
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As a result the density-density correlation function turns out to be 



' E < K " h{"j - h >= -t( J o + A 2 ) (A.22) 



N 2 2 

<»,j> 

Finally, the expectation value we were searching for is 



<H >=UJ2 < K - " J) >= -^Uo + J i) (A.23) 



<i,3> 

Recalling that 



< H x >= -tN^- (A.24) 
da 



altogether the average energy is 



<H>=E = -tN^L - ^(II + if) (A.25) 
Using the identity |^ = Dli(a) the average energy is 

E = -tNDh - ^^(/ 2 + If) (A.26) 
For the sake of simplicity it is convenient to introduce the reduced energy as follows 

The minimization condition = for the reduced energy is equivalent to that of the energy E except 
for the fact that a minimum of E corresponds to a maximum of e. Furthermore, we have 

e (t = 0) = l^E(t = 0) = -^-^ (A.28) 
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B The Asymptotic Formulas for Interacting Fermions in D = 2 Di- 
mensions 



In section 4 we have addressed the issue on the existence of a critical value (t/U) c for which a quantum 
phase transition exists in two spatial dimensions. To achieve this goal in our PV approximation scheme 
we need to evaluate the a — > oo asymptotic behaviour of the basic integrals Iq, I\ in ( [2.29D , ( |2.30 ). This 



was easily done in D = 1 dimensions in section 3 but it becomes much more cumbersome in D = 2 
dimensions and we present in this appendix the main details of this calculation. 

Let us consider the integral Iq in two dimensions as the paradigm of asymptotic evaluation: 



T( \ - r f 2n dkldk2 - m-n 

0laj Jo Jo (2vr) 2 cosh(2a(cosfc 1 + cosA:2)) 1 ' ' 

In the limit a — > oo only the points in the neighbourhood of the region 

cos k\ + cos &?2 = (B.2) 

contribute by a significant amount to the integral ( |B . 1| ) . But these points ( |B.1| ) form precisely the 
Fermi Surface (F.S.) of the system at half-filling. This surface is drawn in Figure 5. It is a square 
inscribed in the Brillouin zone with corners at the points A(±7T, 0), B(0, ±7r). 

In order to isolate this contribution from the Fermi Surface we must distinguish two types of points on 
that curve. 



Bl. Regular Points Contribution to Iq 



We call a regular point of the F.S. to any point such as P r in Figure 5 which is not a corner. Thus 
we have 4 sides of regular points each of which equally contributes to Jo- We restrict our attention to 
one side of the square, say the up-right one. 

The contribution of the regular points come from a rectangle centered at the corresponding F.S. 
portion of infinitesimal width 2A (see Figure 5). We must also isolate the contribution of the corner 
points A and B so that we introduce another infinitesimal quantity e such that the length of the 
rectangle is \^2tt — 2e. 

From Figure 5 it is apparent that both infinitesima are not independent. In fact, 

A = e (B.3) 

They are equal but we keep both for the sake of keeping track clearly where everything comes from. 
Given a regular point on the F.S., we can parametrize any point in the relevant rectangular region of 
Figure 5 by two parameters ko, k as 

k\ = ko + k 

(B.4) 

ki = TT — ko + K 

where ko measures the position on the F.S. and k measures the distance from the F.S., that is , 

k G (e,7r- e) 

(B.5) 

k G (-A, A) 
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Thus, on each fringe the energy function cos k± + cosk 2 contributes like — 2 sin Uq sin n. Substituting 
this into Iq we have, 

f A 1 
dko / dn- 



(2tt) 2 J e J -a cosh (4a sin ko sin k) 

As A C 1 we can approximate sinK ~ k. Changing variables to x = (4asmko)K it is possible to 
perform the a —* oo limit on the boundaries of integration in n to obtain 

1 f 00 1 
dko : — r~ / dx- 



(2tt) 2 J e 4asinkQ J-oo coshx 

These integrals are known and collecting the factor 4 from the same number of fringes we arrive at the 
leading behaviour for Iq D ~ 2 ^ coming from the regular points (R.P.), 

( B - 6 ) 

We arrive at a 1/a behaviour but we must be cautious with it for it comes with a lne term which 
means that it is singular as e — ► 0. This means that equation ( jB,6| ) is not the full story about Iq and 
the contribution from the singular points A and B should contain a singular term which cancels out 
We now check that this is in fact the case. 



B2. Singular Points Contribution (Van Hove Singularity) to Iq 



We have left out for the last moment the contribution from the corner points A(-±ir, 0) and B(0, ±7r). 
This amounts to the integration over the dashed region in Figure 5. Due to the identifications of the 
Brillouin zone boundaries, we can collect the 4 dashed regions into 2 dashed squares centered at A(ir, 0) 
and B(Q,ir) by gluing them together. 

As a matter of fact, both dashed squares contribute the same amount to Iq so that we concentrate only 
on the evaluation of one of them, say A(n, 0). Let us denote by vHa this dashed square around A. We 
can parametrize any point in vHa by K%, k<i as 

k\ = IT + K\ 

(B.7) 

k 2 = k 2 

Thus, on this dashed square the energy function cos k\ + cos k 2 contributes like 

cosfci +cosA;2 ~ 2^ K i _ K ^ 

Notice that this is a tipically van Hove-like singularity behaviour revealing the saddle-point nature 
of the singular points A and B. Inserting flB.8| ) in Iq we get 

' ! dK\dK 2 - 



(2-7r) 2 J JvH A cosh(a(K 2 — 

Now it is convenient to transform to light-cone coordinates k±, 

«i = 7j( K + + K -) 
«2 = 7s( K + - K -) 



(B.9) 
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and thus we arrive at 



^ dn+dn. 



(2tt) 2 Jo Jo cosh(2a:K + K_) 
Moreover, this integral can be evaluated by changing to variables (u, v) as, 

u = K + K_ 

(B.10) 

V = K + /K- 

Eventually, after some algebra we obtain the following contribution from the van Hove point A: 

1 - e --i-, fdu ! n " , (B.ll) 



2vra 2(2tt) 2 Jo cosh(2mt) 

Notice that the lne contribution in ( B.ll|) when multiplied by 2 singular points yields the desired 



cancelation of the singular behaviour obtained in (BJ3) from the regular points 



In all, we are left with a finite contribution that can be extracted from ( |B.11| ). Therefore, the 
leading order behaviour of Io(a) as a — > oo is given by 

(B.X2) 

The remarkable fact about this asymptotic behaviour is the presence of the In a term in addition to 
the 1/a term already found in D = 1 dimensions. This in turn is the distintive characteristic feature 
of the D = 2 dimensional problem. 

The technique outlined above to isolate the leading contribution of Io(a) can also be applied to 
dIo(a)/da and dl\{a)/da (in fact Ii(a) is easier to handle and does not requires that technique). In 
Table 2 we summarize the final results as well as the intermediate results showing explicitly that the 
cancelation of the singular lne terms also occurs for dlo(a) / 'da and dl\{a) / da. 
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C Fokker- Planck Hamiltonians for Fermions 



In this appendix we want to describe in more detail the problems arising when trying to apply the 
FP-method of reference |^(| to fermions. The question we want to answer is whether one can construct 
a hamiltonian for which the ansatz ip{a) (|2.15| ) becomes an exact eigenstate. Following [^] we should 
act with Hq on the state ( 2.15| ) and using the anticommutation relations of the operators nj — 5, cj 



and c\ 



H \if>(a) >= -— E ex P( - f E Tik ~\ ^ T jfc + rest unchanged) |CTW > 

<i,j> <i,k>,k^j <j,k>,kj^i 



= -jE ex P("« E T ^"« E T ifc + | E T im )\CDW> (C.l) 

<i,j> <i,k>,k^=j <j,k>,k^=i <l,m> 

If the operators {Tij} commute among themselves we could take the term ^ J2<i m> Tim to the right of 
the exponential recovering in that fashion the state ip(a). This is precisely what happens in the cases 
described in reference |2(J. Here however the operators do not commute and this complicates matters. 
We can apply the Baker-Campbell-Hausdorf formula to the last exponential in ( |C.1| ). The commutator 
between the sum of the first and second sums with the third gives operators which acting on \CDW) 
vanishes. Hence to order a 2 equation ( p.l[) becomes 



H \i(,( a ) >~ -j E exp(-a E T ^ ~ a E T jk)\^) > (C.2) 

<i,j> <i,k>,k^=j <j,k>,k^=i 



So that one gets up to order a 2 : 



^ FP |^(a)>=0 (C.3) 

with 

Hfp = U E {("* " ~\) + ex P("« E ^ik — a E Tjk)} (C4) 

<i,j> <i,k>,k^=j <j,k>,kj^i 

Expanding the exponential in powers of a we observe that the terms linear in a correspond precisely 
to Ui with the correct relationship between a and t/U ( p.33| ). In fact, from these type of comparison is 
how one derive the a-parameters as functions of the coupling constants. Now, to order a 2 we also find 
the operators appearing in U2, but there are some more operators which acting on the CDW- vacuum 
give nothing. The solution of the PV-equations is by no means unique. 

Higher order parameters in a give some new terms but also others that are already contained in U\ and 
Ui- This means in particular that the relation between t/U and a is rather non-linear, t/U = g(a) = 
(D — i)a + . . .. In order to improve the ansatz one should add to the exponential of the ansatz the 
unwanted operators that appear in the FP-hamiltonian (we mean those operators which do not appear 
in the original ansatz). This program has been successfully implemented for Ising-like models but it 
remains to be seen whether it can also be extended for fermions. 
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Table captions 

Table 1 : Reduced ground state energy e(t/U = 1/2) values for the isotropic case as compared to 
several methods in D=l dimensions. 

Table 2 : Asymptotic leading order behaviour for the basic integrals in D = 2 dimensions showing 
the cancelation of the singular In e terms between the contributions of the regular points on the Fermi 
Surface and the van Hove singularities. 
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Figure captions 



Figures 1 a) - d) : Reduced energy e(a;t/U) in ID eq. fl2.28|) for several values of the coupling 
constant t/U showing the bifurcation phenomena responsible for the existence of a critical point (t/U) c 
in the PV approach. 

Figure 2 : The coupling constant t/U as a function of the PV parameter a for ID (solid line) 
showing the existence of a critical point (t/U) c and 2D (grey line) showing the In a behaviour for large 
a. 

Figure 3 : Staggered magnetization versus the coupling constant {t/U): exact ID result (solid 
line), PV ID result (strong solid line) and PV 2D result (grey line). 

Figures 4 : Reduced energy versus the coupling constant (t/U): exact ID result (black line), PV 
ID result (black-grey solid line) and PV 2D result (grey line). 

Figure 5 : Brillouin zone for 2D fermions at half-filling showing the Fermi Surface formed of 
regular points P r and van Hove singular points A(ir,0) and B(0,ir). A = e are the cut-offs. 
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Method 


Reference 


e(t/U = 1/2) Isotropic case 


Neel 




1 


Free Fermi Sea 


M 


1.6766 


Present PV solution 




1.7232 


Jordan- Wigner transformed UHF solution 


1 


1.7292 


RVB ansatz (M = 3) 


i 


1.7580 


Exact solution 


i 


1.7726 



Table 1: Reduced ground state energy e(t/U = 1/2) values for the isotropic case as compared to several 
methods in D=l dimensions. 



D = 2 Integrals 


a — > oo Behaviour 


R. P. Contribution 


Van Hove Contribution 


!q (a) 


In a 

27TCJ 


-^-lne 

na 


J_ l n e + In" 




In q 


^lne 


^rlne ^ 


/r 2) («) 


(f) 2 






da («) 


In a 
24r> :j 


— ttj~t lne 


1 l , lna 

^ me I ^ 



Table 2: Asymptotic leading order behaviour for the basic integrals in D = 2 dimensions showing the 
cancelation of the singular lne terms between the contributions of the regular points on the Fermi 
Surface and the van Hove singularities. 
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